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Abstract 
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1 Introduction 

The AdS/CFT correspondence [1, 2, 3] provides an important theoretical framework for 
studying a class of strongly coupled quantum fleld theories, and at the same time provides a 
non-perturbative approach to quantum gravity. Traditionally, the correspondence has been 
used to study the dynamics of non-abelian gauge theories, since the prototype examples of the 
correspondence relate string theory/gravity in AdS spacetimes to supersymmetric Yang-Mills 
theories. However, in recent times the correspondence has been extended to study various 
condensed matter systems, especially those in the vicinity of a quantum critical point (see [4] 
and references therein). An example is the possible application to non-relativistic fermions 
at unitarity [5, 6], which is the focus of the current paper. 

Fermionic systems at unitarity, which are realized in experiments with cold atoms, are 
described by a non-relativistic CFT. This fleld theory enjoys a symmetry algebra called 
the Schrodinger algebra [7, 8, 9]. The authors of [5, 6] proposed a class of spacetimes which 
could describe fleld theories with non-relativistic conformal invariance holographically, whose 
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metric is 

ds^ = r' {-2dx+dx--(3^r^dx+y + d^l) + ^. (1.1) 

These have the Schrodinger algebra as an isometry algebra, and we will henceforth refer to 
them as Schr spacetimes. A Galilean CFT in d-spatial dimensions is conjectured to be dual 
to Schrrf+3. 

Subsequently, the Schr5 spacetime was realized in string theory [10, 11, 12].^ It was 
obtained by starting from a known type IIB solution, the AdSs x Xr, spacetime, where 
A's is a Sasaki-Einstein manifold, and applying a solution-generating transformation, either 
the Null Melvin Twist [33, 34] or the TsT transformation [35], to generate the spacetime 
Schrs X A5. The TsT transformation involves a T-duality, followed by a shift, and then a 
second T-duality. This can be used in any spacetime with a U{1) x f/(l) isometry and, in 
general, can be thought of as a Melvin Twist, since the transformation adds a background 
NS-NS B-field, effectively Melvinizing the spacetime. The Null Melvin Twist can be shown to 
be equivalent to the TsT transformation in the special case when one of the U{1) isometries 
is null. 

Since the string-theoretic embedding of the Schrs spacetimes can be achieved by using 
dualities, we can see that the dual field theory is a deformed version of = 4 SYM theory 
(for the special case of = S^). The transformation in [10, 11, 12] breaks the SU{4:) 
symmetry of A/" = 4 down to an SU{3) x f/(l) subgroup, and the field theory is then the 
theory twisted by the R-charge associated with the U{1) isometry. One can also view the 
field theory as the one obtained by decoupling the open string field theory living on D3- 
branes in a Null Melvin spacetime. More generally any A/" = 1 superconformal field theory 
in 3 + 1 spacetime dimensions (corresponding to a dual Sasaki-Einstein compactification) can 
be deformed to a two dimensional (spatial) non-relativistic CFT using the U{1)r symmetry. 

Using the same solution-generating transformation, the authors of [10, 11, 12] constructed 
a two-parameter family of black hole spacetimes with Schrs asymptotics. These black holes 
can be shown to solve the five-dimensional effective equations of motion obtained by a 
Kaluza-Klein reduction of the 10-dimensional type IIB theory on the A5. The simplest five 
dimensional effective action is composed of gravity coupled to a massive vector field and a 
single scalar field [10], which is a truncation of a more general five- dimensional Lagrangian 
involving three scalar fields. The latter, remarkably, is a consistent truncation of type IIB 
supergravity on [11]. 

This black hole solution was used to study the equilibrium thermodynamic properties 
of the field theory and was shown to be dual to the grand canonical ensemble for the dual 
field theory. The thermodynamics is, not surprisingly, consistent with non-relativistic scale 
invariance in two spatial dimensions. In particular, it was found that e = P, where e is the 

^For other recent studies of non-relativistic systems in the context of gauge-gravity duahty see [13, 14, 15, 
16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]. 
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energy density and P is the pressure, as required in non-relativistic CFTs. Furthermore, 
non-equihbrium transport properties of the non-relativistic plasma were also explored in 
[10, 12]. In particular, the shear viscosity rj of the non-relativistic fluid was calculated and 
found to take the universal value rj/s = l/Air typical of strongly interacting field theories 
with gravity duals. 

The present paper is concerned with constructing the gravitational dual of arbitrary fluid 
flows in the non-relativistic CFT. The basic premise is to build on the recent ideas in the 
fluid-gravity correspondence [36] to construct inhomogeneous black hole solutions with Schrs 
asymptotics. The hydrodynamic description of a fluid (either relativistic or non-relativistic) 
is an effective field theory which captures the universal long wavelength physics when the 
system achieves local thermal equilibrium. Building upon previous discussions of the hy- 
drodynamic description of four-dimensional superconformal field theories in the AdS/CFT 
correspondence (see [37] and references therein), in [36] it was argued that starting from the 
most general black hole solution in AdS (a boosted Schwarzschild-AdS black hole) one can 
promote the temperature and velocity fields to local functions of the boundary coordinates. 
It is then possible to solve for the bulk metric order by order in a boundary derivative ex- 
pansion and recover the relativistic Navier-Stokes equations, which encode the conservation 
of the boundary energy-momentum tensor, entirely as a consequence of Einstein's equations. 
Importantly, it is possible to show that the resulting bulk solutions are genuine black holes, 
i.e. that they have a regular event horizon [38]. This correspondence has been extended in 
many directions, including to forced fluids [39], to conformal field theories in various dimen- 
sions [40, 41, 42], to charged fluids [43, 44, 45], to Bjorken flow [46, 47], and incompressible 
non-relativistic fluids [48]. 

We are interested in constructing inhomogeneous black hole solutions with Schrodinger 
asymptotics and extending the fluid-gravity correspondence to relate them to fluid flows in 
the non-relativistic conformal theory. As a first step, we note that a four parameter family 
of black hole solutions with Schrs asymptotics can be obtained by boosting the solutions 
considered in [10, 11, 12]. We can then promote these parameters to functions of the field 
theory coordinates and find the solution order by order in a boundary derivative expansion. 
However, we will show that we can obtain inhomogeneous black hole solutions with the de- 
sired asymptotics in an easier way, by simply applying a TsT transformation to the solutions 
constructed in [36]! We will argue that this in fact captures all the hydrodynamic properties 
of the system in the planar limit. 

To relate these inhomogeneous black hole solutions to fluid flows, we need to calculate 
the boundary stress tensor from the bulk solutions. Because of the slow fall-off of the metric 
perturbation in the asymptotically Schrodinger black holes, the usual technique of obtaining 
the stress tensor by functionally differentiating the action with respect to the boundary 
metric cannot be straightforwardly applied to these cases (see Appendix A for a discussion 
of this issue). However, [11] proposed that the stress tensor of the asymptotically AdS 
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space before the TsT transformation can be re-interpreted as the stress tensor complex of 
the non-relativistic theory. We will adopt this approach. We therefore describe the general 
reduction of a relativistic stress tensor on the light cone to obtain a non-relativistic stress 
tensor complex. The structure of the relativistic conformal stress tensor implies that for any 
non-relativistic conformal theory obtained in this way, the thermal conductivity k of the 
non-relativistic fluid is 

e + P , , 

K = 2r]^- 1.2 
pT 

where p is the mass density, or, even more succinctly, 

Pr = 1 (1.3) 

where Pr is the Prandtl number. 

The outline of this paper is as follows: we will begin in §2 by discussing how the rel- 
ativistic fluid equations are reduced on the light-cone to the non-relativistic Navier-Stokes 
equations. This will allow us to explore the general properties of the non-relativistic stress- 
tensor complex. We will review some aspects of the Schrs solutions in §3 and then describe 
how to construct inhomogeneous black hole solutions dual to arbitrary fluid flow and give 
the dual stress tensor to first order in derivatives in 54. We end with a discussion in 55. 



2 Light-cone reduction of relativistic fluids 

Consider a relativistic fluid in Minkowski space in + 2 spacetime dimensions; we will use 
light-cone coordinates {x^,x~,x.} and take the metric to be 

dSflat = Vt^u dx^ dx" = -2 dx~ + dx^. (2.1) 

Suppose we view this fluid in the light-cone frame and evolve it in light-cone time x"*". 
Then, for fixed light-cone momentum P„, we obtain a system in d + 1 dimensions with non- 
relativistic invariance. This is of course familiar from the discrete light-cone quantization 
(DLCQ) of quantum field theories. In fact, one of the models for studying non-relativistic 
conformal field theories holographically, suggested in refs. [13, 14], was that one could con- 
sider pure AdS, with the relativistic conformal symmetry broken to Galilean symmetry 
simply by compactification of the x~ coordinate, which singles out a preferred light-cone di- 
rection. Note that in this case we are not only compactifying the light-cone direction in the 
boundary where gravity is non-dynamical (and the metric flat, (2.1)), we are also required 
to compactify the coordinate in the bulk AdS spacetime. This involves introducing closed 
null curves in the geometry and the validity of supergravity becomes questionable [11]. We 
will return to a different gravitational background, viz. (1.1), where x" does not need to be 
compactified to achieve Galilean symmetry. Note however that it is still useful to take x~ 
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compact, so that the momentum P_ is integer quantized, since P_ is interpreted as particle 
number in the dual theory. 

Relativistic hydrodynamics in + 2 dimensions is formulated in terms of pressure (or, 
equivalently, the temperature) and the four velocity u^, subject to the condition that ri^^, = 
— 1. This gives d + 2 degrees of freedom. At the same time, non- relativistic hydrodynamics 
in d spatial, one temporal dimensions can be formulated in terms of the mass density p, the 
pressure P, and the spatial velocities v\ also giving d + 2 degrees of freedom. 

We would like to find a mapping between the degrees of freedom of the ((i+2)-dimensional 
theory to the degrees of freedom of the d + 1 dimensional theory such that the relativistic 
hydrodynamic equations imply the non-relativistic hydrodynamic equations. We would also 
like to find how the thermodynamic quantities of the two formulations are related. Finally, 
we plan to use the map to find out the thermal conductivity of the non-relativistic theory. 
We will first begin with ideal hydrodynamics and then discuss dissipative terms. 

2.1 Ideal fluids 

The relativistic hydrodynamics equations are just the conservation of energy and momentum 

V^T^"" = 0. (2.2) 

An ideal relativistic fluid has a stress tensor given by^ 

T^" = (erei + Prci) + Prci V^" , (2.3) 

where the energy density erei is related to the pressure Prei by a thermodynamic equation of 
state. Equations (2.2) and (2.3) define a system of d + 2 equations for the d + 2 unknowns. 

Non-relativistic ideal hydrodynamics is described by the continuity equation, 

dtp + d,{pv') =0, (2.4) 
together with the equation of momentum conservation, (here i = 1, . . . ,d) 

dt{pv') + djli'^ = 0, rf^' = pv' + , (2.5) 
and the equation of energy conservation, 

dt(^e + ^p v'^ + d, jl = 0, jl = ^{e + P)v'v\ (2.6) 

where f ^ = f * f 

^We use the subscript "rel" for quantities in relativistic hydrodynamics and indicate quantities in non- 
relativistic hydrodynamics without subscripts. 
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Consider the relativistic equations (2.2) on the hght-cone. We will consider only solutions 
to the relativistic equations that do not depend on x~] that is, all derivatives d- vanish. The 
coordinate x"*" corresponds to the non-relativistic time t. The equations of energy-momentum 
conservation are, 

(9+T++ + diT+' = , (9+T+^ + djT^ = , d+T+- + diT-' = 0, (2.7) 

which reduce to the non-relativistic equations under the following identification: identify T"^"*" 
with the mass density, with the mass flux (which is equal to the momentum density), 
T*-' with the stress tensor, T"* with the energy density, and T~* with the energy flux, 

T++ = p, T+' = pv\ T^ = W^^ 

T+- = e + ]^pv\ T-''=jI. (2.8) 

It is now easy to convince oneself based on (2.8) that the precise mapping between 
relativistic and non-relativistic hydrodynamic variables is 

= \ - — - — , u' = u^v\ 
\ 2 e + P ' 

Prcl = P, e,ei = 2£ + P. (2.9) 

The component of the relativistic velocity u~ can be determined using the normalization 
condition u^u'^ = —1 to be 

""=K^^ 

While the analysis has been for a general relativistic fluid with an equation of state 
erei(-PrGi); we wiU soou focus on conformal fluids. Conformal invariance requires that the 
stress tensor for the relativistic theory be traceless, = 0, which gives us the equation of 
state erei = (rf + 1) Pfci- In the non-relativistic theory we can once again use the conformal 
invariance to learn that 2e = dP. 



2.2 Viscous fluids 

We now wish to extend our mapping of relativistic hydrodynamics into non-relativistic hy- 
drodynamics to first order in derivatives on both sides. The ideal stress energy tensor (2.3) 
can be supplemented with dissipative terms, which can be expanded systematically in terms 
of derivatives of the velocity field and pressure. Specifically, we have 

T'^" = (e,el + Prel) W + T/^'^P.el + V^^^ (2.11) 

where n^'^ incorporates all the dissipative contributions. For first order viscous hydrodynam- 
ics we have 

7^'^'^ = -2r/relr'^^ (2.12) 



6 



where 

T^" = \ P'^ {^cUf, + ^f,u^ - ^ i]^f, V^u^ (2.13) 

is the shear tensor and we have introduced the spatial projector P^''^ = rj^'^ + u'^ . 

We will use the zeroth-order equations of motion to simplify the viscosity term. By 
using zeroth-order equations, we make an error of second order in derivatives, which can be 
neglected. Namely, we use the ideal hydrodynamic equations in the following form, 

W^u" + ^^^p = 0, where = P'^'Va , (2.14) 

rcl 

to rewrite the stress-energy tensor as 

T"" = (€„, + F„,,)i,''i,'- + F„,<'' 

V d+1 erel + Prel J 

On the non-relativistic side, we use the ideal hydrodynamic equations in the form 

dtp + d,{pv')=0, 

dtv' + div' + -diP = 0, 
P 

dte + di{ev')+PdjV^ = 0. (2.16) 
The first-order contributions to the (spatial) stress tensor and the energy flux are 

ji^ = [s + P + \p v'^) + ri a'^ - K diT, (2.17) 

where k is the thermal conductivity and T is the temperature. 

By using (2.15) and (2.17), we now establish the mapping between relativistic and non- 
relativistic viscous hydrodynamics. First, we find r^^ = 0, and therefore 

T++ = (e,ei + Prci) (m+)1 (2.18) 
The identification T~^^ = p implies then that 



n+ = J-^— , (2.19) 

fcrol ~r J rol 



unchanged from the ideal hydrodynamic level (2.9). 



Next, we find 
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On the other hand, we still want to map = pv^. This means that there is now a 
correction to the relation between and f*: 



= U' 



p \ 2 (ercl + Prcl) 



(2.21) 



For T*-' , after some algebra, we find 

=pv' + Prel - Vrel U+ {d,V, + d,Vi - ^ 5'^ S^V^^ , (2.22) 

which implies that the pressures on the two sides still coincide, 

Prcl = P, (2.23) 
and the relationship between the viscosities is 

Vrel = 4- (2-24) 

Note that our identifications automatically give a first-order correction cTjj in the non- 
relativistic theory with the correct tensor structure. That is, the trace-free relativistic shear 
tensor gives a trace free spatial stress tensor in the non-relativistic theory. 

Regarding the other components of the stress tensor, after some calculations involving 
many cancellations, one discovers that is 

T+- = ^(erei-Prei) + ^pw', (2.25) 

which means that the relationship between relativistic and non-relativistic energy densities 
remain unchanged, 

ere\ = 2e + P. (2.26) 

Finally, for T~* we find 

/ 1 \ (9-7/+ ?/ + 

T-'=[e + P+-p - r/,ei u+a'^ v, + r/,ei S'' - Vrd — 5'' d,P. (2.27) 



Thus we have to require that 



+ 



Vrel 7^ - Vrcl — d,P = - K d,T . (2.28) 

In order to see that the left hand side is indeed proportional to SjT, we need to use the 
mapping (2.19) and the equation of state for a non-relativistic theory. 

Focusing specifically now on conformally invariant fiuids, using (2.19) and e = ^ P, we 
find 

diu+ M+ le + P / p{d+4)/(d+2)\ 

Vrcl T^TT? - ^7^1 — diP = -r^rci \ / <9i In . (2.29) 

[u+y p V 2p VP/ 



Recalling that the equation of state of the holographic non-relativistic liquid is [25] 



{d+2)/2 



P = <^[-) . (2.30) 

the argument of the logarithm in (2.29) is up to a constant. Therefore, the left hand 
side of (2.28) is indeed proportional to 9jT, and one reads out the value for the thermal 
conductivity: 

e + P 

K = 27]^-. (2.31) 
pT 

Let us now compute the Prandtl number. The Prandtl number is defined as the ratio of 
the kinematic viscosity u and the thermal diffusivity 

Pr = -, (2.32) 

X 

where 

^ = ^ , X = — , (2.33) 

P pCp 

where Cp is the specific heat at constant pressure. We note the definition of the heat capacity 
at constant volume: 

where H = E + PV is the enthalpy. Write H = wV = w N/n. We then find 

Cp = N—(-) = :r{l^] (2-35) 



dT \n J p in? \ dT J 



p 



where we have used the fact that w = {d/2 + 1) P and is fixed at fixed P. At fixed P, fi ^ T^, 
and n = dP/dT ~ and dn/dT = —2n/T. Therefore 



Cp = (2.36) 



2wN 
Tn 

and Cp = Cp/ M (M being the total mass) is equal to 2w / pT. Thus we find: 



Pr = — i = 1. (2.37) 

pi K 

Note that this result is valid for any non-relativistic conformal fluid obtained from the DLCQ 
of a relativistic conformal fluid. 



3 Thermal description of non-relativistic CFTs 

Non-relativistic conformal field theories with Schrodinger symmetry in d spatial dimensions 
are dual to Schr^+s spacetime [5, 6]. The Schrs spacetimes can be realized in string theory as 
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the near- horizon geometry of D3-branes probing a Null Melvin universe [10].^ Furthermore, 
this background can be obtained as a solution to a 5-dimensional Lagrangian which is a 
consistent truncation of IIB supergravity [11], involving gravity coupled to a massive vector 
field and three scalars: 

5buik = tttV / d'x^ (R + V{<P,)-5{d<P,r-^-^{d<P2r-l{dhf 
lovrGs J \ 2 2 

V{(t)i) = 24e-'^^-^'^2 -4e-^'^i-^'^^ -8e~^°'^^ 

In fact, all known solutions are solutions to a slightly simpler theory outlined in [10], where 
the three scalars are linearly related as 

{01,02,03} = {-|-^,1}0. (3.2) 

This allows one to consider the five-dimensional effective action 

S = J d'x^ (^R - ^-{d,<^){d^<^) - h-'^''F^,F>^^ - AA^A^ - , (3.3) 

where 

\/(0) = 4e2'^/3(g2<^_4)_ (34) 

The field theory on the D-branes in the appropriate decoupling limit is A/" = 4 Super 
Yang-Mills plus a non-commutative deformation, giving rise to a version of the dipole field 
theories discussed in [33, 49]. Consider the fields in A/" = 4 SYM (with gauge group SU{N)) 
which transform under the global symmetry S0{A,2) x S0{6). Picking a U{1) C S0{6) we 
deform the field theory by replacing ordinary products appearing in the Lagrangian by a 
star-product [11] 

f^9 = e^'^'-'''-'-^'h9 (3.5) 

where P_ is the momentum along the null x~ direction and R is the f/(l) charge. This prod- 
uct is a hybrid between the usual non- commutative star-product which only involves spatial 
momenta and the /3-deformation of A/" = 4 which involves only the R-charges. Of importance 
later will be the fact that this deformed field theory inherits some of the properties from 
A/" = 4 SYM. In the large limit we in fact expect that the planar sector of the deformed 
theory to be identical to that of A/" = 4 SYM [11].^ 

•^The near-horizon limit in this case needs to be taken keeping in mind that we want to reahze the scahng 
symmetry consistent with the Schrodinger algebra. We thank James Lucietti for a useful discussion about 
this issue. 

''As discussed in §1 these statements extend to generic J\f = 1 superconformal field theories which we deform 
by a U{1)r symmetry. 
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The Lagrangian (3.3) has a black hole solution which was obtained by the Null Melvin 
Twist [10, 12] or TsT solution-generating transformation [11] in 10 dimensions. The black 
hole geometry is given as 

- r'^ fir) 



ds% = r'^k{r) 3 



fir) 



4/3 



2 



+ \ idx-)^ - [1 + fir)] dx+ dx' 



+ fc(r)^ + (3.6) 

with the massive vector and scalars taking the form 

kir) V 2 
1 



(3.7) 



where /(r) and /c(r) are 



/W = l-4' MO = H-^, (3.8) 



^^ = W7^2■ (3.10) 



with 7^ = /3^r^. Note that in these light-cone coordinates, the solution asymptotically 
approaches the vacuum solution (1.1) at large r, and also reduces to (1.1) when we set 
r+ = 0. It has been argued in [10] that the black hole spacetime (3.6) corresponds in the 
field theory to a grand canonical ensemble at temperature 

T = ^ (3.9) 
and a chemical potential for particle number 

1 

Note that the Schrodinger algebra involves a conserved charge associated with particle num- 
ber, which geometrically is realized via the Killing field These results were obtained with 
minor differences in the derivation in [11, 12]. 

Using a Euclidean action calculation, the authors of [10] derived the conserved charges 
of the black hole. This calculation required a careful analysis of the boundary counter- 
terms, since the metric (3.6) has rather complicated asymptotics. In Appendix A we discuss 
the calculation of the conserved charges in a Hamiltonian formulation, supplementing the 
analysis of [10]. The conserved charges associated with the Killing symmetries d- and d+ 
in (3.6) translate in the field theory to particle number N and total energy E. To obtain 
finite values we assume that the ^"-direction is compactified with period Ax~ . The results 
obtained in [10] are 
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and 

where we have given the results in terms of the physical temperature and chemical potential, 
and also in terms of the parameters in the bulk solution to emphasize the physical interpre- 
tation of the parameters 7 and r_|_. Furthermore, the pressure is given in the grand canonical 
ensemble directly in terms of the Gibbs potential Q(T, V^): 

PV^-Q(T,,,V)^^^V^.-, (3.13) 

leading thus to an equation of state 

PV = E =^ P = e. (3.14) 
In addition, we note for future reference that the entropy of the black hole is given by 

lo G5 



4 Hydrodynamic description of non-relativistic CFTs 

The black hole solution (3.6) corresponds to an equilibrium configuration of the non-relativistic 
field theory. We would like to study departures from equilibrium in the continuum limit by 
considering local patches of the field theory in local equilibrium. This is the hydrodynamic 
limit, where one has local variations of energy density and particle number; these local 
domains evolve according to the laws of hydrodynamics which were described in §2. 

To obtain a gravitational description of the fluid dynamical regime, we need to patch to- 
gether local domains of equilibrated fluid - the precise manner in which this can be achieved 
was outlined in [36]. The authors presented an algorithmic procedure to construct gravita- 
tional solutions starting from the equilibrium black hole solution. The idea is to consider the 
most general stationary solution for the equilibrium - for the case of relativistic supercon- 
formal theories in c? -|- 2 dimensions, this is given by a boosted Schwarzschild-AdS^+s black 
hole, which is specified by c? -|- 1 parameters: a horizon size r_|_ and a unit normalized velocity 
field Ufj_. Choosing Eddington-Finkelstein like ingoing coordinates which are regular on the 
horizon, the metric takes the form: 

ds^ = 2u^ dx^dr - /(r) u^u^ dx^'dx'' + P^^ dx^dx" (4.1) 

where P^u = Vfiu + is the spatial projector introduced earlier and /(r) = 1 — {r+/rY~^^. 

To study the hydrodynamic configurations one promotes the parameters r+ and to 
functions of the boundary coordinates x'^. Then one recursively solves for the bulk metric 
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order by order in a boundary derivative approximation. This procedure was systematically 
carried out to second order in [36]. To ensure regularity of the bulk solution, and in particular 
to guarantee the existence of a regular future event horizon, it was necessary to adapt 
coordinates wherein lines of constant boundary x'^ corresponded to ingoing null geodesies 
in the bulk. This coordinate chart was utilized in identifying how the locally equilibrated 
boundary domains evolve in the bulk radial direction - the fluid in such a domain was evolved 
along a tube in the bulk centered about the ingoing null geodesic. The size of these tubular 
domains is set by the local energy density and they provide the bulk analog of patching 
together pieces of equilibrated fluid. 



4.1 Gravitational dual of a non-relativistic fluid: Direct construction 

We now turn to the analogous calculation for the non-relativistic CFT discussed in §3. We 
should first write the metric (3.6) in a form that is regular through the future horizon. To 
do so, perform a coordinate transformation 



X ' 



x"*" + I3p{r) 



X 



(4.2) 



with p{r) chosen such that the metric is regular on the horizon r = r_|_. A convenient choice 

5 

P'(^) = ^TTT ' (4-3) 



IS 



leading to a metric 

ds^ = r'^ k{r) 



1 - f{r) 2 ft ; 



{dx^ 
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{dx~ 



[1 + /(r)] dx+ dx- 



2/5r2 dx-^ + 2(3dx'' 



dr - k -i 0^ dr^ + k{r) 3 r^dx^ 



(4.4) 



This gives a two-parameter family of solutions of the 5d effective Lagrangian (3.3). We 
can construct a four parameter family of solutions trivially by performing a boost via the 
coordinate transformation 



X + VX , X X +V-X+-t' X . 



(4.5) 



Since the background metric (1.1) has Galilean invariance, this boost does not change the 
asymptotic form of the metric. 

This boosted black hole, characterized by the parameters {r_|_, /5, v}, is the starting point 
for a systematic hydrodynamic analysis. Following [38] we should promote the parameters 



^If one wanted further dj. to be null then one could instead choose p'{r) = PTI^' 



13 



{r+,/3,v} to fields depending on (x"'~,x) and consistently solve the equations of motion 
order by order in derivatives in the and x directions. This process, while straightforward, 
is rendered cumbersome by the presence of additional fields in the action (3.1). We will 
therefore resort to a trick to recover the dual geometry. 

4.2 Inhomogeneous black holes via TsT transform 

The solution we want to obtain is characterized by four parameters depending on the coordi- 
nates {x~^, x). The general inhomogeneous black hole of [36] also depends on four parameters, 
r_|_ and u^, and if we specialize to solutions independent of x~ , they also depend on the coor- 
dinates (a;"*", x). For such solutions, we can then apply the same TsT transformation used to 
obtain the black hole solution (3.6) to obtain inhomogeneous black holes with Schrodinger 
asymptotics. Since these solutions will by construction reduce to (3.6) when the parameters 
are constants, we argue that they are precisely the required inhomogeneous solutions. Below, 
we discuss the action of the TsT transformation on the solution of [36]. 

The TsT transformation is a solution-generating technique in string theory wherein one 
uses a twisted T-duality to add NS-NS flux to a given solution [35]. Consider a solution to 
type IIB supergravity of the form Ai^ x where is a Sasaki-Einstein space, which we 
view as a U{1) fibration over a base B, and Ai^ is a constant negatively curved spacetime, 
which we will take to be asymptotically AdSs. The simplest case is when A5 = and 
the base is thus a CP^. We assume that Ai^ admits a Killing vector field d-, so that 
A^s X A's has an isometry group U{l)x- x f/(l)^. Viewed as an eight dimensional solution 
upon Kaluza-Klein reduction, we have an SL{2, R) symmetry group, which can be used to 
generate new solutions. This is the TsT transformation. 

We start from the solutions dual to arbitrary fluid flow for hydrodynamics of relativistic 
superconformal theories of [36], 

ds^ = dX^ dX^ = -2 Uf, S{r, x) dx^ dr + x,,u{r, x) dx" dx" , (4.6) 

where X^ = {x^, r}. The metric dual to viscous fluid dynamics for the relativistic conformal 
fluids is given by (4.6) with S{r, x) = 1 and 

2 2 

Xfiu = {Pfiu - f{r, x) Uy) H F(r, x) t^^ + - r V xu^ - r u^V x [uyU^) , (4.7) 

f+ "J 

where P^y is the spatial projector defined after (2.13), f{r,x) is the function /(r) appearing 
in (3.8) with r+ — > r+(x), and 

F{r,x) = ^ 

It is important to note that and r_|_ are no longer parameters, but regarded as functions of 
x^. To apply the TsT transformation, we need to choose a direction x~ , and assume r+, 
are independent of x~. 
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In 



(r + r+(x)) (r + r_|_(x)^ 



2 arctan(r/r+(x)) + vr 



(4. 



The general map between the parameters of a relativistic fluid and the parameters of 
a non-relativistic fluid to first order in derivatives was written in (2.19), (2.21) and (2.23). 
When we write the relativistic fluid in terms of the gravitational dual spacetime, we can 
rewrite this mapping as a mapping between the functions r+, characterizing the asymp- 
totically AdS spacetime and the functions r+,/3, appropriate in the non-relativistic case. 
The function r+ is the same in both descriptions, and the map from to (3, is 



(4.9) 



The full 10 dimensional metric is a direct sum of the metric qab on Ai^ given in (4.6) 
and a Sasaki-Einstein space As, 

dsl = QAB dX^ dX^ + (# + Af + ds'^{B), 

F(5) = 4 (VoKA^s) + hYo\{B) A (# + A)) . (4.10) 

Under the TsT transformation the metric (4.10) gets mapped to a new solution of Type IIB 
super gravity [11]^ 

ds\ = {gAB dX^ dX^ - e^"^ qa^ Qb^ dX^dX^ + e^^ {dtp + Af + ds\B)) , 

F(5) = 4 (Vol(A^5)+Vol(i3) A/i(# + A)), 
5(2) = e'^'^h'^gB^dX^ Aidtfj + A), 

= l + h^g_^ . (4.11) 

This solution can be Kaluza-Klein reduced back to five dimensions to give a metric (restrict- 
ing to situations where the norm of the Reeb vector is fixed to = 1), 



dsl = e-y {gAB-e^'^gA-gs-) dX^ dX"" , (4.12) 



supported by 



,2 



A = e^'f'gA-dX' 
1 



1+g— 



(4.13) 



As before this five dimensional solution solves the equations of motion arising from (3.1) 
with the scalars being related as in (3.2). 



^In [11], the TsT transformation involved an arbitrary parameter a; however, this can be absorbed into 
a redefinition of the coordinate x" by a boost in the plane. We have found it more transparent to fix 
the parameter in the TsT transformation and keep instead the velocity /3 in (4.9) as the free parameter 
corresponding to the choice of x" . 
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Restricting to configurations which have d- as an isometry, after a TsT transformation 
on the metric (4.6) we get a new metric of the form 





2u^S dx^ dr + xab 



1 + 



Xa Xb 



dX^ dX 



B 



A 



e^'t'XAdX^, 



1 



(4.14) 




with 



Xa = Xa~ - 5 5^ . 



(4.15) 



The TsT transform converts the asymptotically AdSs spacetime (4.6) to an asymptotically 
Schrs spacetime, which depends on the r_|_, (3, Vi defined in (4.9) which are arbitrary functions 
of (x"*", x). This provides the required inhomogeneous generalization of the black hole solution 
(3.6). We will not write the result of applying the transformation more explicitly, as it is 
quite complicated, and its construction is a straightforward exercise. Also note that since 
the internal manifold plays a minimal role in our construction, it is easy to verify that the 
solution (4.14) also solves the consistent truncation action (3.1) with the scalars still related 
according to (3.2). 

4.3 Properties of black holes dual to non-relativistic fluids 

We will now discuss some of the physical properties of the solution (4.14). The solutions 
(4.14) are the most general long- wavelength regular solutions dual to configurations of the 
dual non-relativisitic conformal field theory and are valid to leading order in the boundary 
derivative expansion. These geometries solve the field equations arising from (3.3) provided 
the boundary stress tensor complex satisfies the non-relativistic Navier-Stokes equations. 
This follows from the fact that (4.6) are the most general regular solutions dual to the 
relativistic field theory and that the regularity properties of the black hole solutions are 
unaffected by the TsT transformation. The regularity of the solutions (4.6), in particular, 
the fact that they have a regular event horizon, was demonstrated in [38]. This result required 
that the variations in the boundary directions parameterized by are slow. In employing 
the TsT transformation, all we required was that d- be a Killing vector; so there is no 
variation of the relativistic fluid in the light cone direction and one can of course take the 
variations in all the other directions to be appropriately slow. To explicitly demonstrate that 
the solution (4.14) has a regular event horizon, one can follow the perturbative construction 
of [38]. From this analysis it is easy to infer that for viscous non-relativistic fluids, the 
location of the horizon will remain at r = r_|_(x'^,x). 

In analyzing the regularity of the geometries constructed in [36], it was important to 
work in a well behaved coordinate chart. As explained there and subsequently elaborated in 
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[38, 42] the solutions can be thought of as being tubewise approximated by a homogeneous 
black hole solution, with the tubes being domains in the bulk centered around radially 
ingoing null geodesies with width set by the scale of variation in the boundary directions x^. 
In the non-relativistic case (4.14), the metric is not a priori written in coordinates adapted 
to the radially ingoing geodesies. To see this note that the gauge choice employed in [36, 38] 
was to set Qrfj, oc m^.^ On the other hand, in the non-relativistic case we have non-trivial grr 
arising after the TsT transformation coming from the non- vanishing Xr- This can of course 
be removed by a coordinate transformation, and then the issue of regularity boils down to 
the analysis presented in [38]. 

To complete the fluid-gravity correspondence, we need to identify the stress tensor com- 
plex associated with the inhomogeneous solutions constructed above. The boundary field 
theory dual to the Schrs background was argued to be a non- commutative deformation of 
A/" = 4 SYM. This fact allows us to argue that the planar sector of the deformed theory with 
non-relativistic invariance is in fact identical to the parent J\f = A theory [11], and the direct 
computations of the thermodynamics in [10] are consistent with this picture. We therefore 
argue that, as advocated in [11], we can identify the non-relativistic stress tensor complex 
corresponding to the geometry after the TsT transformation with the one we have prior to 
the transformation. That is, the dual stress tensor complex is obtained simply by performing 
the DLCQ reduction described in §2 for the relativistic stress tensor corresponding to the 
geometry (4.6) before the TsT transformation. 

More explicitly, the relativistic stress tensor corresponding to (4.6) (to first order in 
derivatives) is of the form (2.11) with® 

erel = 3Prel, Prel = , Vrel = ■ (4-16) 

167rG5 167rG5 

Thus, using the identifications in section 2, the non-relativistic stress tensor complex dual 
to the geometry (4.14) will be of the form (2.17) with^ 

The boundary stress tensor complex for the solution (4.14) comprises a spatial stress 
tensor Iljj, particle density p, an energy flux j*, momentum flux pv^ and energy e. While 
we have argued that in the large planar limit these quantities can be derived by reducing 



^In fact, as discussed in [42], one can simplify the metric further, by demanding that the radially ingoing 
null geodesies are afhnely parameterized by r. This is equivalent to setting S(r,x) = 1 in (4.6). 
^Note that ^q^q^ gives the effective central charge of the dual field theory. For the case of deformed A/" = 4 

SYM to the non-rclativistic dipolc theory, this evaluates to where N is the rank of the gauge group. 

^Note that the stress tensor complex obtained from the relativistic stress tensor in §2 is a local density in 
the x~ direction. From the non-relativistic field theory point of view, it is more natural to multiply by Ax~ 
to obtain an object which is a local density only in the spatial directions. 
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the relativistic stress tensor on the hght-cone (using thus the properties of the TsT trans- 
formation), it should in principle also be possible to obtain these by direct computation in 
the geometry (4.14). It turns out to be easy to use the counter-term construction proposed 
in [10] to extract the spatial stress tensor without trouble. However, it is not clear how 
to calculate the energy e and the particle density p in an analogous fashion. In [10], the 
energy e and the particle density p for the system in thermal equilibrium were obtained by a 
Euclidean action calculation. In Appendix A, we discuss the calculation of these quantities 
using canonical methods. 

5 Discussion 

We have discussed the hydrodynamic limit of d spatial dimensional non-relativistic confor- 
mal field theories and their dual gravitational solutions, which are inhomogeneous black 
holes with Schr^+s asymptotics. We employed the fact that starting with a relativistic 
hydrodynamical system one can obtain non-relativistic fluid dynamical equations upon an 
appropriate light-cone reduction. Relativistic hydrodynamics m. d + 2 spacetime dimensions 
descends to non-relativistic hydrodynamics in d spatial dimensions. We demonstrated this 
both for ideal fluids and also for dissipative fluids where we restricted attention to first order 
in the gradient expansion. The latter allowed us to recover the heat conductivity of the 
non-relativistic system in terms of the shear viscosity and state parameters of the parent 
relativistic fiuid. In particular, for non-relativistic CFTs we have shown that the Prandtl 
number is unity. ^° Since the first order dissipative coefficients for the non-relativistic CFTs 
arise from the shear viscosity of the relativistic system, it is not surprising that the rates of 
momentum diffusion and thermal conduction are correlated. 

Our construction of the geometries dual to the non-relativistic fluids uses as its starting 
point the asymptotically AdS spacetimes dual to relativistic hydrodynamics. Starting with 
the solutions constructed in [36] we were able to construct inhomogeneous black holes with 
Schr asymptotics using the TsT transformation described in [11]. The TsT transform deforms 
the boundary fleld theory to a non-local quantum fleld theory through the introduction of 
a star-product (3.5) - we are therefore looking at the hydrodynamic limit of these deformed 
superconformal fleld theories. 

Given these black hole solutions we can in principle try to extract the non-relativistic 
stress tensor complex of the dual fleld theory using an appropriate boundary stress tensor 
construction. Instead of doing so, we adopted the philosophy advocated in [11] - the hydro- 
dynamics of the non-relativistic theory is given by the light-cone reduction of the undeformed 
superconformal theory in the planar limit. This can be justifled by recalling that for super- 
conformal fleld theories with gauge group of rank A^, the star-product deformation resulting 

^''For comparision, Pr(watcr) ~ 7, Pr(mcrcury) ~ 10""^, and Pr(air) w 0.7. 
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from the TsT transformation leaves the planar sector of the theory unchanged [35]. This 
allows us to efficiently extract the stress tensor for the non-relativistic viscous fluid. One 
should in principle explicitly compute the stress-tensor complex using the counter-term sub- 
traction scheme discussed in [10]. It turns out to be easy to check that one can indeed extract 
the spatial stress tensor. Furthermore, by rewriting the action in a Hamiltonian formulation 
we have been able to extract the energy and particle density. But we have encountered some 
difficulties in proving that the Hamiltonian generates the time-translation symmetry of the 
non-relativistic field theory. This is an interesting open problem which we hope to return to 
in the future. It is also worth remarking that to compute equilibrium thermodynamics we 
can take a simpler route proposed in [25] - use a background subtraction scheme where the 
reference spacetime is not the vacuum Schr^ spacetime, but corresponds to a state with zero 
temperature at finite particle density. 

Our discussion of the light-cone reduction of the relativistic equations was restricted to 
ideal and viscous fluids, i.e., up to and including the first order in the gradient expansion. 
Fluid dynamics in general can be viewed as an effective field theory with an infinite number 
of irrelevant terms obtained as usual in a derivative expansion. It would be interesting to 
understand the light-cone reduction of the second order relativistic conformal hydrodynamic 
stress tensor constructed in [36, 50] (see [41, 42] for the result in arbitrary dimensions). 
In particular, the various relaxation times encountered at second order should descend to 
interesting transport coefficients for the non-relativistic theory. It must however be men- 
tioned that unlike the relativistic case, where causality issues [50] require us to consider 
non-linear hydrodyanamics, one is not similarly forced to consider higher-order terms in the 
non-relativistic setting. It would also be interesting to carry out a systematic analysis of 
non-relativistic conformal hydrodynamics by looking at the constraints on the allowed ten- 
sor structures coming from the Schrodinger symmetry, paralleling the relativistic analyses of 
[50, 51]. 

The fluid dynamics we have discussed here has been restricted to conformal fluids with 
Schrodinger symmetry. One main feature of such fluids is that they are in general compress- 
ible; this follows from the fact that the energy density is related to the pressure 2e = dP 
through the equation of state (which in turn follows from scale invariance). To make contact 
with the usual studies of incompressible Navier-Stokes equations we need to ensure that we 
can decouple the fluctuations in the density. This can be achieved by looking at low fre- 
quency modes which do not excite the propagating sound mode in a hydrodynamic system. 
In fact, this limit was discussed recently in the context of the fluid-gravity correspondence in 
[48, 52] where the authors showed that starting from a parent relativistic conformal fluid dy- 
namical system one can recover incompressible Navier-Stokes equations in a suitable scaling 
limit. Curiously the limit procedure reveals an interesting structure in the fluid equations - 
they are scale invariant under a new scaling symmetry. This symmetry is different from the 

^^In the notation used in [10] this corresponds to the hmit r+ — > and /3 ~y oo with 7 = /3r5|_ held fixed. 
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Schrodinger symmetry enjoyed by the fluids under consideration in this paper. We arrived 
at our hydrodjTiamic description in d spatial dimensions via a light-cone reduction of a + 2 
dimensional relativistic theory. In [48, 52] the authors derive non-relativistic incompressible 
hydrodynamics m d + 1 spatial dimensions by a suitable limit of the relativistic theory. It 
would be interesting to the take the incompressible limit of the non-relativistic conformal 
fluid described here to compare with their results. 
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A Hamiltonian calculation in the bulk 

Since the field theory dual of the asymptotically Schr geometries is non-relativistic, we might 
not expect the construction of a spacetime stress tensor for asymptotically AdS geometries 
to have a natural counterpart for these solutions. We might instead expect the map from 
bulk to boundary to involve some kind of Hamiltonian framework, which treats spatial and 
time directions separately. In this appendix, we describe the construction of a Hamiltonian 
for the bulk geometry from the action constructed in [10]. This is a generalization of the 
calculation of [53], using the known form of the action to determine the correct boundary 
terms for the Hamiltonian framework. 

The idea of the calculation in [53] is to choose a foliation of the spacetime by constant 
time slices, and rewrite the covariant action in the form S = j dt[pq — H], keeping careful 
track of the boundary terms. This will then give us a form for the Hamiltonian including 
boundary terms at infinity. For simplicity we start from the truncated version of the action 
used in [10], 

S = J d'xV^ (^R - ^{d,<p){d>^<j>) - ^-e-"f/'F,,F^'' -4A,A^ - V{<p)^ (A.l) 

+— ^ / d^^v^ {2K -Q + A^A^' + cA^Ay + b{A^A^y + {2c-Ab + 3)(l)^) 
lovrGs J 

for some arbitrary constants c, d. To simplify the formulae below, we will work with c = 6 = 0, 
but it is straightforward to carry out the analysis in the general case. We take the constant 

^^The analysis can be easily extended to the consistent truncation Lagrangian of [11] modulo the cost of 
some cumbersome expressions. We find that we have an eight parameter set of boundary terms which lead 
to the finite on-shell action reported in [10]. 
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time slices to be the surfaces 11+ of constant x~^. The coordinates on S+ will be denoted 
collectively as x'^ and the induced metric is k^^}^ The lapse and shift A^, A^'^ are defined by 
considering a vector u'^ defined on S4., such that V^x"*" = 1, and decomposing this vector 
into the normal to S_l and the shift, 



(A.2) 



We also have a constant r surface which is our cut-off on the spacetime, with coordinates 
and induced metric h^^^. To top off the list of surfaces we have the constant time surfaces on 
the cut-off surfaces which we will denote as and we reserve for the coordinates and 



a^iy for the induced metric on S^^. 



The analysis of the gravitational part of the action (A.2) is very similar to the calculation 
in [53], and leads to 



A^ 



K 



(A.3) 



where p^jy is the gravitational conjugate momentum, p^j/ = ^Jl^{K^l, — Kg^^)^ and r'^ is 
the unit normal to and n = det(K^jy). The contributions to the bulk Hamiltonian and 
Hamiltonian density constraints are 



Ha 



167r Gr, 



3 ^/K 



Here is the covariant derivative associated with the spatial metric k^u- Only the boundary 
term will contribute to the on-shell Hamiltonian. We note that this can be rewritten as 



/r 
rf'C v^Y^(2(%'-6Ar)-2Ar'^p; 



(A.5) 



For the matter part of the action, 
1 , — / 4 



Sr, 



] 

4 



3 ^ 



performing the space-time split and introducing the momenta 

1 



1671 G. N 



(A.7) 
(A.8) 
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We will not distinguish between the indices used for different hypersurfaces. 
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where dot denotes differentiation wrt x~^, we can rewrite this action as 

2v^ 



+ p\A^ + 



fi,m 



+ / dx^ 



1 

IQttG.N 



N 



(A.9) 



We need to integrate the term (c)^yl_|_) by parts; will then be a non-dynamical field, 
and we can eliminate it using its equation of motion, 



A. = mA,, -271 G.N V, 



This gives us the matter action in its final form. 



Pa 



dx~^ I d'^x 

2+ 



+ 



dx^ 
N 



dX^/^ 



+ p\Ax-Nn^-mn^,^ 
1 



(A.IO) 



(A.li: 



Pa 



IottGsA' Jk 



167r Gf 



where 



Stt G^^/H 



V, 



+ 



5„2 



P; + 



Stt a 



5 89i/3 M „v 



(A.12) 



167r 



and 



(A.13) 



Thus, the matter contribution to the on-shell Hamiltonian is 

1 



On-shell 



Pa 



IQtx G.N 



{A+-N'^A^Y-A+^r^ + 



N 



K 



IQ-kG, 



A, - N'^Ar. oc A^ 



0, 



(A.14) 
(A.15) 



so we can drop the first term, which will not contribute to the value or first variation of the 
Hamiltonian, and write 



d^C 



(A.16) 



which is closer in form to the gravitational part. 
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For our spacetime (3.6), we have (writing for brevity i{r) = \/JJr) k{ry^^) 

d 



9x+ 



= 7 d r(l + /(r)) d ^ 



which leads to 



AT = -1^3 £(r), 

7 



,l + f{r) d 



and r'^ being the unit-normal to is just given as 

F(r) d 
(r) dr 



r' = "4^l- (A.19) 



Direct computation gives 



Vfi = 7r, 

(3)^ = 

i{r) ' 

p_r N 7^ 



,/k SttG, r^f{r)k{ry 

The gravitational part of the Hamiltonian then involves 



/ ^ 72 21 7^ _ 



(A.20) 



We see that the divergence cancels. The divergence will cancel against a contribution 
from the matter part, which we turn to next. 



On the 3-boundary, 



72 



SO 



A^A^ + 3<p'^^--^. (A.23) 

12 



We also evaluate 
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to obtain 

/ 72 21 7^ 



5 



^^^^--L^^Cv/^^ 5i,--i,). (A.25) 



As a result, the total Hamiltonian is 

H"' = tttV / d'C^N% = --^ / rf=^C<, (A.26) 

in agreement with our expectations: this matches the total energy obtained by Euclidean 
methods in [10]. 

It is easy to also extract the particle number density from this calculation; we simply 
consider shifting the vector ^ + a^, where a^d^t = 0, which shifts A^^ N^^ + a^. 
This will change 

and hence redefines the Hamiltonian by 



(A.28) 

We should interpret this new Hamiltonian associated with + a'^ as a combination of the 
energy and the momentum in the direction specified by a^. We see immediately that if 
a'^ points in one of the spacelike directions, the change in the Hamiltonian vanishes, which 
agrees with our expectation that the spatial momentum densities vanish in the state we are 
considering. If we take a'^ = a5'^__, then the gravity term is 

P^r^ _JL_ 1' rfjr) ^ f_ 
^ SttG^ f{r)k{r) £(r) SyrGgr^' ^ ' ^ 



The matter terms give 



e"^*^/^ aN 1 o A^ -y^ 

a^F.^r'^^--^!., T,-^ KA.)^ - ^ ^, (A.30) 



so only the gravitational term contributes to the change in the Hamiltonian, which is 

if- ^ - / d\ v^iV-^ 1^ = ^^- I d'C (A.31) 

consistent with the value (-P-) = 7^^Ax~/47rG'5 obtained in the Euclidean approach. 

Thus, this approach gives us a definition of the Hamiltonian whose on-shell value appears 
correct. We should ask if the Hamiltonian so defined is the generator of the asymptotic time- 
translation symmetry on this class of spacetimes. To be the generator of the symmetry, our 
Hamiltonian should satisfy [54] 

5if = y d\ [5p^, E^^ + E"^ + bk^'' A^, + 6^ , (A.32) 
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with no surface terms, where ifj denotes the collection of matter fields. Now, our Hamiltonian 
constructed from the covariant action should satisfy this property by construction, since 
(A. 32) is just the requirement for Hamilton's principle 6 j dt [pq — H) = to have solutions, 
and the boundary terms in our action S = j dt [pq — H) were chosen precisely so as to ensure 
that SS = on-shell. However, in [10], the action was only shown to satisfy 6S = for a 
restricted set of boundary conditions on the variations, relating the leading order variations 
of different fields. It would therefore be useful to check directly that the Hamiltonian we 
have proposed satisfies (A. 32). We have encountered difficulties in performing this check; we 
leave their resolution to future work. 
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